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SUMMARY 


This  Report  determines  the  accuracy  with  which  the  Seasat  1  altimeter 
height  measurement  can  be  used  to  derive  the  spacecraft  CG  to  mean  sea 
level  distance.  To  calculate  this  accuracy  a  description  of  the  mean  shape 
of  a  reflected  pulse  must  be  used  which  incorporates  a  realistic  descrip¬ 
tion  of  a  wavefield’s  waveheight  probability  distribution  in  its  derivation. 
A  quantitative  description  of  such  a  reflected  pulse  is  developed  here. 

This  description  is  then  used  in  a  simulation  procedure  to  calculate  values 
of  the  rms  errors  in  the  tracking  and  wave  biases,  hence  the  rms  height 
error.  This  figure  of  accuracy,  ie  rms  height  error  is  derived  for  various 
sea  state  conditions.  _ 
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1  INTRODUCTION 

The  Seasat  1  satellite  launched  in  May  1978  was  the  first  satellite  devoted  to 
observation  of  the  World's  oceans.  To  perform  this  task  one  passive  and  three  active 
microwave  instruments  were  carried  onboard  this  satellite.  These  instruments  consisted 
of  a  synthetic  aperture  radar  used  for  producing  two-dimensional  images  of  the  ocean 
surface,  a  scanning  multichannel  microwave  radiometer  from  which  sea  surface  temperature, 
wind  speed  and  atmospheric  water  vapour  concentration  could  be  derived  and  a  scat terome ter 
and  altimeter.  The  scatterometer  was  basically  a  radar  system  used  to  measure  the  sea 
surface  wind  vector  while  the  altimeter  measured  the  rms  waveheight  of  those  waves  lying 
in  a  region  of  sea  directly  beneath  the  satellite  and  the  satellite's  height  above  the 
sea  surface. 

With  this  instrument  complement,  the  sea  state  variables  of  rms  waveheight  and 
dominant  wavelength  and  direction  were  measured  along  with  the  meteorologically  significant 
variable  of  the  sea  surface  wind  vector  in  all  weathers. 

One  further  instrument  was  carried  -  a  visual  and  infra-red  radiometer  -  which  gave 
information  about  the  prevailing  weather  conditions  in  the  field  of  view  of  the  four 
microwave  instruments. 

This  Report  is  concerned  with  the  altimeter  instrument.  As  mentioned  above  this 
instrument  measured  the  rms  waveheight  and  satellite  altitude.  To  derive  the  altitude 
the  transit  time  of  a  very  short  duration  pulse  was  measured.  From  this  transit  time 
measurement  the  distance  between  the  mean  sea  level  and  the  satellite's  CG  (centre  of 
gravity)  can  be  calculated.  Combining  this  distance  with  the  geocentric  distance  of  the 
satellite  the  distance  corresponding  to  the  sum  of  the  geoid  and  instantaneous  sea  surface 
height  can  be  found.  This  quantity  is  of  considerable  interest  to  oceanographers  and 
geodesists.  By  determining  this  sum  both  spatially  over  the  oceans  and  in  time,  oceano¬ 
graphic  features  such  as  mesoscale  eddies,  ocean  tides  and  currents  can  be  measured. 
Included  in  the  latter  are  the  currents  due  to  the  large  scale  ocean  circulation  patterns. 
This  circulation  is  responsible  for  heat  transport  from  the  tropics  to  the  poles  -  conse¬ 
quently  it  has  an  effect  on  the  climate.  Currents  can  be  measured  by  an  altimeter  since 
a  current  introduces  a  slope  into  a  sea  surface.  This  slope  can  be  detected  by  measuring 
the  change  in  the  above  sum  along  the  spacecraft's  ground  track.  Unfortunately  the 
gradient  of  this  slope  is  very  small  -  typically  only  a  few  centimetres  of  height  varia¬ 
tion  occurs  over  a  horizontal  distance  of  a  hundred  kilometres.  Consequently  the  accuracy 
to  which  the  satellite’s  altitude  can  be  found  must  be  as  high  as  possible  if  such  a 
small  slope  is  to  be  detected. 

To  derive  the  satellite's  altitude  a  series  of  corrections  must  be  applied.  These 
corrections  are  known  as  the  instrumental  and  geophysical  corrections.  The  instrumental 
corrections  basically  translate  the  satellite  reference  to  the  satellite  CG  while  the 
geophysical  corrections  correct  for  the  propagation  delay  introduced  by  the  ionosphere 
and  troposphere. 

In  this  Report  it  is  assumed  that  the  instrument  and  geophysical  corrections 
necessary  to  correct  the  measurement  of  altitude  made  by  the  Seasat  I  altimeter  have  been 
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applied'.  The  instrument  corrections  are  accurately  known  and  the  uncertainty  in  the  geo¬ 
physical  corrections  is  less  than  a  few  centimetres,  provided  that  radiometer  data  is 
available  to  determine  the  magnitude  of  the  wet  component  part  of  the  tropospheric  correc¬ 
tion  to  the  measured  altitude.  The  error  sources  to  be  investigated  here  are  due  to 
tracking,  discriminator  and  wave-bias  noise.  They  arise  due  to  the  presence  of  noise  on  a 
pulse  received  after  reflection  by  a  sea  surface  and  to  the  error  this  noise  introduces 
into  the  values  of  the  rms  waveheight  and  skewness  variables.  Their  magnitudes  are 
obtained  by  fitting  a  theoretical  description  of  the  mean  shape  of  a  received  pulse 
envelope  to  that  obtained  in  practice. 

To  determine  the  accuracy  of  the  Seasat  I  altimeter  altitude  measurements  under 
varying  sea  state  conditions  (a  range  of  rms  waveheight  and  skewness  values) ,  a  realistic 
description  of  the  mean  shape  of  a  received  pulse  envelope  (hereafter  known  as  a  mean 
waveform)  must  be  developed.  The  theoretical  basis  of  such  a  description  is  discussed  in 
the  first  section  of  this  Report  along  with  a  review  of  how  the  Seasat  1  altimeter 
measured  the  spacecraft's  altitude  by  tracking  to  a  specific  point  on  a  mean  waveform's 
leading  edge.  The  remainder  of  the  Report  is  concerned  with  using  this  theoretical 
description  of  the  mean  waveform  to  determine  the  wave  and  tracking-bias  errors  and  in 
conjunction  with  tracking  error,  the  accuracy  of  an  altitude  measurement  made  by  this 
altimeter. 


2  ALTIMETER  DESIGN  AND  OPERATION 
2. 1  General 

Up  to  the  end  of  1980  only  three  spacecraft  had  included  an  altimeter  within  their 
instrument  complement.  These  spacecraft  were  Skylab  -  the  first  American  manned  space 
laboratory,  GEOS  3  a  satellite  used  for  geophysical  research,  and  Seasat  I,  the  first 
satellite  devoted  primarily  to  oceanographic  remote  sensing.  Of  these  three  altimeters, 
the  one  carried  by  Seasat  1  measured  the  altitude  with  the  lowest  level  of  instrumental 
range  noise  (10  cm  as  opposed  to  around  50  cm  for  the  altimeter  on  GEOS  3  and  approximately 
I  m  for  the  Skylab  altimeter) .  To  determine  the  accuracy  in  the  measurement  of  the 
distance  between  the  spacecraft  CG  and  mean  sea  level,  derived  from  a  Seasat  1  altimeter, 
height  measurement,  an  understanding  of  how  this  altimeter  measured  the  satellite's 
altitude  is  needed.  The  design  and  operation  of  this  altimeter  will  therefore  be 
described. 


2.2  The  Seasat  I  altimeter 

The  design  of  this  altimeter  is  illustrated  schematically  in  Fig  I  and  its  para¬ 
meters  are  given  in  Table  I . 

To  operate  successfully  this  system  has  first  to  acquire  the  backscattared  signal 

in  the  acquisition  mode  and  then  remain  locked  onto  it  in  the  tracking  mode.  In  the 

acquisition  mode  an  initial  measurement  of  the  transit  time  of  a  pulse  is  made.  Owing 

*3 

to  the  low  orbital  eccentricity  (e  <  10  )  of  Seasat  1  the  altitude  only  varied  between 
±40  ka  centred  on  800  km.  Consequently  the  transit  time  falls  within  the  range  5.06  to 
5.6  as.  To  measure  this  tiae  a  monochromatic  pulse  of  duration  3.2  pa,  frequency 
13.5  GHz  is  transmitted.  At  the  time  of  triggering  this  pulse  a  range  counter  within  the 


82059 


5 


SACU  (synchroniser,  acquisition  and  control  unit)  counts  down  for  a  pre-selected  time 
interval  before  triggering  the  production  of  a  CM  LO  (local  oscillator)  pulse  of  frequency 
13.0  GHz  and  duration  3.2  us.  After  first  stage  mixing  the  resulting  signal  is  detected 
and  passed  to  a  threshold  detector  within  the  SACU.  If  this  threshold  is  exceeded  the 
value  of  the  pre-selected  time  interval  is  used  to  form  a  'height  word',  ie  a  word  holding 
the  best  estimate  of  the  transit  time  at  that  point.  This  word  is  then  given  to  the  ATU 
(adaptive  tracking  unit)  and  the  tracking  mode  initiated.  However,  if  the  threshold  is 
not  exceeded,  further  monochromatic  pulses  are  transmitted  with  the  range  counter  counting 
down  through  a  series  of  pre-determined  time  intervals  until  the  threshold  is  exceeded  and 
acquisition  attained. 

The  tracking  mode  differs  fundamentally  from  the  acquisition  mode  owing  to  its  use 
of  the  pulse  compression  technique.  The  principle  advantage  of  this  technique  lies  in 
its  reduction  of  the  peak  transmitted  power  for  the  attainment  of  a  specific  SNR  (signal- 
to-noise  ratio)  by  eliminating  the  need  for  very  short  (~ns)  duration  pulses.  As  knowledge 
of  this  technique  is  essential  for  understanding  the  operation  of  the  altimeter  in  the 
tracking  mode  the  pulse  compression  technique  will  now  be  discussed. 

2.3  The  pulse  compression  technique 

In  this  technique  a  chirped  pulse,  ie  a  pulse  whose  carrier  frequency  is  swept 
linearly  from  13.5  to  13.82  GHz  in  3.2  us  (in  the  case  of  Seasat  I)  is  employed.  Consider¬ 
ing  only  reflection  by  one  point  on  a  surface  distance  H  from  the  altimeter  (point  B 
in  Fig  2),  then  2  H/C  seconds  after  transmission  of  the  pulse  leading  edge,  the  leading 
edge  of  the  reflected  pulse  is  received.  On  mixing  this  received  pulse  with  a  LO  FH 
(frequency  modulated)  pulse  whose  frequency  sweep  rate  equals  that  of  the  transmitted 
pulse  but  whose  starting  frequency  is  13  GHz,  a  constant  difference  frequency  F  (equal 
to  500  MHz  in  the  case  of  Seasat  1)  for  a  specific  height  H  will  be  observed.  If  H 
varies  F  will  vary.  Employing  this  procedure,  all  points  at  the  same  range  will 
produce  mixer  outputs  at  the  same  frequency.  Consequently  integrating  the  power  from  the 
mixer  at  a  particular  frequency  will  give  the  total  power  backscattered  and  received  by 
the  altimeter's  antenna  from  those  points  on  a  sea  surface  at  the  same  range. 

Onboard  the  spacecraft  the  received  signal  is  sampled  by  a  set  of  60  contiguous 

filters  which  comprise  the  DFB  (digital  filter  bank).  The  swept  frequency  rate  is 

100  MHz/us  and  the  time  resolution  of  the  system  is  3.125  ns.  The  filter*  bandwidth  must 
-3 

be  <100  x  3.125  *  10  MHz  {ie  312.5  kHz)  so  as  not  to  degrade  the  system  time 
resolution.  By  integrating  the  DFB  outputs  over  contiguous  frequency  intervals,  each  of 
bandwidth  312.5  kHz,  the  response  from  successive  concentric  annuli  on  a  flat  sea  surface 
is  obtained.  Plotting  the  received  intensity  as  a  function  of  time  gives  the  received 
pulse  power  envelope.  The  resulting  shape  of  this  received  pulse  would  be  precisely  the 
same  if  instead  of  a  3.2  us  chirped  pulse  a  pulse  of  3.125  ns  duration  were  transmitted 
and  the  received  pulse  were  sampled  with  a  filter  of  width  3.125  ns.  The  production  of 
such  a  short  duration  pulse  is  a  great  problem  due  to  the  peak  power  required  to  provide 
a  useful  SNR  ( signal- to-noise  ratio).  As  a  chirped  pulse  width  is  1024  times  greater 


*  Each  filter  is  also  known  as  a  gate 
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than  a  pulse  of  duration  3.J25  ns  its  peak  power  can  be  1/1024  of  that  of  a  3.125  ns  wide 
pulse  to  achieve  the  same  SNR. 

2.4  Tracking  mode  operation 

When  used  in  the  tracking  mode  the  transit  time  is  determined  by  the  ATU.  In 
practice  50  successive  DFB  outputs  are  stored  and  averaged  in  the  ATU .  The  ATU  then  com¬ 
bines  these  averaged  gate  values  in  one  of  five  ways.  The  procedure  is  illustrated  in 
Fig  3.  Combining  the  gates  together  by  this  means  forms  a  tracking  gate.  Five  tracking 
gates  can  be  produced  by  the  ATU,  each  tracking  gate  consisting  of  an  early,  middle  and 
late  gate.  For  each  tracking  gate  the  widths  of  the  early,  middle  and  late  gates  are  the 
same,  the  output  of  these  gates  being  obtained  by  sunning  the  outputs  of  adjacent  gates 
in  the  DFB  which  lie  in  the  early,  middle  or  late  gate  time  range.  The  number  of  DFB 
gates  lying  in  this  range  depends  on  which  of  the  five  tracking  gates  is  being  synthesised. 
This  summing  process  is  carried  out  by  an  algorithm  run  on  a  small  micro— computer  in  the 
ATU.  A  variety  of  tracking  gates  were  used  so  that  the  magnitude  of  the  attitude  error 
correction  to  the  altitude  measurement  was  minimised. 

For  this  system  to  measure  the  transit  time  to  the  mean  sea  surface  the  intensity 

of  that  point  on  a  received  waveform's  leading  edge  corresponding  to  the  return  from  the 

2 

mean  sea  surface  must  be  known.  According  to  the  model  of  Miller  and  Brown  *  (where  the 
waveheight  pdf  is  normally  distributed)  which  describes  the  mean  waveform,  this  intensity 
equals  the  power  level  of  a  mean  waveform  curve  at  the  origin  of  the  time  axis.  This 
intensity  shows  a  small  variation  with  a  this  being  0.443  for  a  normalised  mean 
waveform  .  (In  practice  the  60  DFB  gate  outputs  sunned  to  a  constant  value,  this  being 
achieved  by  the  use  of  an  AGO  system.  Consequently  all  calculations  involving  a  mean 
waveform  assume  it  to  be  normalised.) 

The  time  offset  between  this  level  (whose  normalised  power  level  was  0.443)  and  the 
level  corresponding  to  the  mean  sea  surface  was  calculated  using  the  Miller  and  Brown 
model,  taking  Er  -  0°  and  the  results  tabulated  as  a  function  of  o w  ,  being  stored  in 
the  altimeter's  micro-computer .  By  measuring  the  gradient  of  the  mean  waveform's  leading 
edge,  by  subtracting  the  late  gate  output  from  the  early  gate,  o^  could  be  determined 
and  thus  the  necessary  time  offset  derived  by  interpolation.  Consequently  the  transit 
time  could  be  corrected. 

To  track  to  any  particular  point  on  a  mean  waveform's  leading  edge  a  tracking  law 
of  the  form: 

I(t)  -  M(t)  -  oL(t  «■  W) 

is  used.  Here  M  and  L  are  the  power  outputs  of  the  middle  and  late  gates,  t  the 
offset  of  the  centre  of  the  siiddle  gate  from  the  desired  tracking  location,  and  W  is 
the  width  of  the  middle  and  late  gates. 


*  In  this  model  all  time  measurements  are  swde  with  respect  to  the  tisw  of  reception  of 
the  mean  sea  level  contribution  to  a  waveform. 
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In  practice  the  device  is  tracking  successfully  when  I(t)  «  0  .  By  employing 
variable  gate  widths  the  stability  and  sensitivity  of  the  system  is  maximised  .  The  con¬ 
dition  I(t)  >0  is  attained  as  follows. 

Having  obtained  a  timing  offset  the  ATU  now  tries  to  vary  the  sequencing  time, 
ie  time  between  triggering  the  DDL  (digital  delay  line)  for  the  production  of  the  trans¬ 
mitted  and  L0  chirped  pulses  so  as  to  centre  the  point  on  a  mean  waveform's  leading  edge 
whose  normalised  intensity  is  0.443  in  the  DFB.  The  value  of  this  sequencing  time  is 
expressed  in  the  17  MSB  (most-significant-bits)  of  the  25  bit  'height  word*.  The  value 
of  these  17  MSB  is  used  as  the  countdown  time  by  the  DDG  (digital  delay  generator)  which 
triggers  the  DDL.  The  value  attributable  to  the  8  LSB  (least-significant-bits)  of  this 
height  word  is  used  to  vary  the  centre  frequency  of  the  individual  gates  of  the  DFB  by 
the  same  amount  so  as  to  decrease  the  timing  offset.  This  height  word  is  formed  for  every 
received  pulse,  its  value  being  derived  from  the  previous  height  word  added  to  a  value 
from  a  feedback  loop. 

3 

The  feedback  loop  is  also  part  of  the  ATU.  It  produces  a  timing  offset  by  taking 
the  average  offset  for  the  previous  50  pulse  grouping  and  adds  a  certain  fraction  of  the 
timing  offset  to  a  rate  accumulator  and  another  fraction  to  a  position  accumulator.  The 
sum  obtained  from  these  accumulators  at  millisecond  intervals  then  provides  the  input  to 
form  a  new  height  word.  By  constantly  centring  the  DFB  on  the  pulse  leading  edge, 
ie  fine  tuning,  and  varying  the  sequencing  time,  ie  coarse  tuning,  the  height  word  is 
continually  updated  and  the  device  tracks  smoothly.  In  this  way  the  tracking  condition 
I(t)  -0  is  obtained. 

This  updating  is  needed  because  of  an  altitude  rate  or  acceleration.  When 
travelling  over  the  sea  the  altitude  rate  varied  smoothly,  this  variation  being  due 
mainly  to  orbital  eccentricity.  Over  land  this  rate  was  often  non-linear  due  to  the 
presence  of  underlying  rough  terrain. 

The  DDG  is  contained  within  the  SACU  and  as  mentioned  above  triggers  the  production 
of  a  chirped  pulse.  The  chirped  pulses  produced  for  transmission  and  mixing  differ  in 
their  power  level  and  frequency.  Initially  a  chirped  pulse  centred  on  250  MHz  whose 
frequency  is  decreasing  linearly  in  time  is  produced  by  applying  a  12.5  ns  pulse  of  centre 
frequency  250  MHz  to  a  surface  acoustic  wave  device.  The  resulting  chirped  pulse  has  a 
bandwidth  of  80  MHz  and  a  sweep  rate  of  3.2  us.  The  bandwidth  of  this  pulse  is  then 
increased  by  a  factor  of  four.  For  transmission  the  centre  frequency  of  this  chirped 
pulse  is  up-converted  to  3375  MHz  and  then  multiplied  by  4,  resulting  in  a  13.5  GHz 
chirped  pulse.  For  the  L0  chirped  pulse  the  associated  frequencies  are  3250  MHz  and 
13  GHz. 

In  pra-eice  when  lock,  ie  I(t)  ■  0  has  been  obtained  the  reflected  and  LO  chirped 
pulses  ar  combined  in  the  first  mixer.  An  output  centred  at  500  MHz  is  produced. 

Further  mixing  then  occurs,  with  a  GW  oscillation  of  frequency  500  MHz.  Frost  this  mixer 
two  signals  are  obtained,  the  in-phase  I  and  quadrature  Q  outputs  (Q  is  the  same  as  I 
except  for  a  phase  change  of  90°) .  Both  these  signals  are  then  sampled  by  a  HSWS  (high 
speed  wave  sampler)  which  undertakes  analogue-to-digital  conversion  at  a  rate  of  20  MHz, 

64  samples  being  taken  at  equal  intervals  during  3.2  us. 
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To  complete  the  system,  filtering  by  the  DFB  takes  place  in  those  intervals  between 
successive  pulse  receptions.  For  an  irradiated  sea  surface,  contributions  to  the  received 
waveform  arise  from  various  ranges.  As  each  range  is  associated  with  a  particular 
frequency  the  received  waveform  is  synthesised  from  components  covering  a  range  of 
frequencies.  To  determine  a  waveform  the  power  within  each  gate  must  be  determined  at 
each  sampling  time.  Integration  over  a  period  of  3.2  ps  for  all  the  gates  will  give  the 
waveform.  This  waveform  is  then  taken  by  the  ATU  and  used  in  the  tracking  process. 

Apart  from  performing  this  tracking  operation  the  ATU  also  provides  information  on 
the  rms  wave height,  the  spacecraft's  altitude,  the  AGO  gain  value  and  the  waveform  shape. 
The  average  of  approximately  100  consecutive  measurements  of  rms  waveheight,  AGO  value 
and  altitude  were  telemetered  to  ground  along  with  the  average  of  100  consecutive  waveform 
shapes . 

3  QUALITATIVE  DESCRIPTION  OF  A  MEAN  WAVEFORM 

3. 1  General 

This  description  will  be  developed  as  follows.  Firstly  the  flat  sea  impulse 
response  will  be  derived.  This  response  is  a  description  of  the  waveform  received  when 
a  pulse  of  zero  time  duration  has  been  scattered  by  a  flat  sea  surface.  The  advantage 
of  employing  an  impulse  response  to  characterise  a  surface  lies  in  its  use  in  determining 
the  mean  waveform  for  rough  sea  surfaces  which  are  illuminated  by  pulses  of  finite  dura¬ 
tion.  For  example  if  Ps(t)  describes  the  transmitted  pulse  shape  as  a  function  of  time 
and  W(t)  describes  the  distribution  of  surface  elevations  from  a  mean  level,  then  the 
mean  waveform  P  is  obtained  by  evaluating  the  expression  Ir(t)  *  Ps(t)  *  W(x) 
where  Ir(x)  is  the  flat  sea  impulse  response  and  *  denotes  a  convolution  process. 

To  derive  the  flat  sea  response  the  assumption  is  made  that  the  mean  power  levels 
received  from  elements  of  a  scattering  surface  each  having  unit  area  are  incoherent,  so 
the  powers  sum  linearly.  This  assumption  is  examined  in  section  3.2.  To  derive  the  mean 
power  received  per  unit  area  of  backscattering  surface,  Es  the  scattered  far  field 
electric  vector  must  be  determined.  E8  itself  is  produced  by  a  surface  current  density 
J  that  has  been  induced  in  the  sea  surface  by  H  (the  magnetic  field  vector)  of  an 
incident  electromagnetic  (EM)  wave. 

3.2  Determination  of  the  scattered  far  field  electric  vector  and  backscatter  coefficient 

In  determining  Eg  an  important  criterion  is  whether  the  ratio  R  ,  between  the 
characteristic  scale  size  and  EM  wavelength,  quantified  by  Barrick^  as  2wh(x,y)/\£  is 
greater  than  unity.  Here  h(x,y)  is  the  height  of  point  (x,y)  on  the  sea  surface  above 
mean  sea  level  and  X.  is  the  wavelength  of  the  impinging  radiation.  This  ratio  is 
important  as  the  method  of  deriving  Eg  differs  for  the  two  cases.  The  two  physical 
situations  are  illustrated  in  Fig  4.  The  top  diagram  shows  the  case  R  >  1  .  Here  the 
surface  is  smoothly  varying  on  a  scale  comparable  with  the  wavelength.  If  visible  light 
for  exanple  is  incident  in  the  direction  shown  and  the  backscattered  component  is  viewed 
then  the  bulk  of  the  reflected  radiation  will  be  seen  to  come  from  a  series  of  small 
surface  elements,  the  reflection  being  specular,  is  mirror-like.  This  situation  is  well 
illustrated  by  the  glinting  phenomena  seen  when  sunlight  is  reflected  by  a  water  surface. 
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In  the  region  of  these  elements  -  so  called  'specular  points'  the  surface  electric  current 
density  induced  by  the  impinging  EM  field  must  be  producing  a  scattered  electric  field 
component,  the  phase  of  which  is  only  varying  slowly  across  such  a  region  resulting  in 
constructive  interference.  For  this  to  be  so  the  surface  would  be  expected  to  be  normal 
to  the  impinging  wave.  Further  from  these  regions,  due  to  the  tilt  of  the  water  surface, 
the  phase  variation  between  adjacent  scattered  components  is  much  greater,  so  they  tend  to 
cancel.  This  picture  is  the  same  for  microwave  radiation  with  wavelengths  of  the  order  of 
centimetres  incident  on  a  smoothly  varying  surface  of  scale  size  greater  than  tens  of 
centimetres. 

The  scattered  electric  field  component  Es  can  be  solved  by  the  use  of  an  integral 
expression^  derived  directly  from  Maxwell's  equations.  This  takes  the  form: 

(I  -  kgkg)  •  J(cr')  exp(ikr  -  iks  •  r')ds' 

where  u>  is  the  angular  frequency,  k  the  wave  number  of  the  impinging  EM  wave  and  ks 
the  wave  number  vector  of  the  scattered  wave.  1  is  the  unit  dyad  and  J  the  surface 
current  density.  The  integral  is  evaluated  over  the  actual  scattering  surface,  all 
primed  symbols  being  measured  on  this  surface. 

To  determine  Eg  a  description  of  J  on  the  surface  must  be  provided.  At  a 
frequency  of  13.5  GHz  seawater  is  a  good  conductor.  Consequently  the  displacement 
current  D  is  much  smaller  than  the  conduction  current,  so  from  Maxwell's  equations 
linking  H,  J  and  3D/3t,  V  A  H  ■  J  to  a  good  approximation.  Adopting  a  coordinate 
system  defined  by  the  three  orthogonal  unit  vectors  a,  g  and  ^  where  a  and  g  lie 
on  the  plane  of  the  mean  surface  and  j  is  anti-parallel  to  k  and  taking  H  »  aH  then: 


T  »  3H  3H 
-  "  2  Tw  -  1  3ft 


The  current  into  the  surface  is  of  no  interest  as  this  cannot  produce  a  detectable 
scattered  Hs  component  as  Hs  is  orthogonal  to  J  ,  consequently  only  the  component 
of  J  parallel  to  the  surface  is  of  interest. 

In  practice  the  so-called  optical  distribution  is  taken  for  J  where  J  ■  24  A  . 
Here  n  is  the  unit  vector  normal  to  the  surface  and  Hi  the  magnetic  field  vector  of 
the  incident  EM  wave.  In  practice  this  form  of  J  is  somewhat  unrealistic.  For  example 
if  a  sea  surface  as  illustrated  in  Fig  5  is  illuminated  as  shown,  regions  of  the  surface 
labelled  'a'  are  not  illuminated,  thus  J  here  must  be  zero.  One  problem  is  what 
happens  to  the  distribution  of  J  at  the  edges  of  the  illuminated  regions.  The  abrupt 
discontinuities  as  implied  by  this  optical  distribution  will  not  occur  as  the  conduction 
current  will  be  expected  to  flow  out  into  the  non- illumine ted  region,  becoming  dissipated 
by  Joule  heating  and  radiation.  Secondly  if  deviations  from  the  mean  surface  occur  and 
the  radi-  on  impinges  at  low  angles  of  incidence  multiple  scattering,  is  scattering  of 
a  art  \  of  the  incident  wave  that  has  already  bean  scattered,  will  take  place  with 
increasing  regularity  (Fig  5). 
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The  two  points  mentioned  ebove  will  complicete  the  distribution  of  J  •  No  descrip— 
tion  of  J  has  yet  been  developed  which  takes  these  two  effects  into  account. 

Normally  altimeters  illuminate  the  sea  surface  such  that  the  angle  of  incidence  to 
the  mean  surface  is  90°.  Consequently  little  obscuration  and  multiple  scattering  is 
likely  to  occur  -  except  in  the  case  of  highly  agitated  seas  where  the  waves  are  breaking. 
Such  sea  conditions  occur  infrequently  as  they  are  associated  with  severe  storms. 

The  picture  is  further  complicated  as  a  sea  surface  is  composed  of  features  ranging 
from  scale  sizes  of  centimetres  to  hundreds  of  metres.  Such  scale  sizes  can  cause  R  to 
be  less  than  and  greater  than  unity  at  different  points  on  the  surface.  Having  both 
roughness  scales  on  a  sea  surface  produces  an  example  of  what  is  known  in  the  literature 
as  a  composite  surface.  The  presence  of  small  scale  structure  (R  <  1  case.  Fig  4) 
produces  diffraction  at  the  surface  -  this  complicates  the  theoretical  treatment  as  to 
evaluate  the  scattering  integral  the  induced  surface  current  density  and  the  shape  of 
the  surface  must  be  known  precisely. 

The  problem  of  scattering  by  composite  surfaces  has  been  examined  by  a  number  of 

9 

authors  (see  Ref  8  and  references  therein)  with  the  mathematical  technique  of  Burrows 
used  by  Brown1^*  being  the  most  satisfactory  to  date  as  it  derives  o*  for  all  angles  of 
incidence  rather  than  just  for  large  and  small  angles  of  incidence  as  other  techniques 
have.  The  essential  result  of  Brown's  work  is  that  for  angles  of  incidence  <5°  the  con¬ 
tribution  to  a  due  to  diffraction  is  under  2%  of  that  of  large  scale  (R  >  1)  scatter. 
This  diffraction  contribution  to  specular  scatter  will  be  ignored  here. 

As  mentioned  previously  the  bulk  of  the  scattered  radiation  comes  from  a  large 
number  of  physically  separate  elements  -  the  specular  points.  To  determine  the  value  of 
E8  at  distance  r  from  a  scattering  surface  one  needs  to  sum  over  the  contributions  of 
each  of  the  components.  As  shown  by  Kodis2  the  contribution  Ejs  to  Es  by  one  specular 
point  for  an  incident  electric  field  of  unit  amplitude  is  given  by: 

Ejs  "  7^  "~'J~  '  elkrf  exP  ik(3  *  £>dydx 

Sj  *  £z  { 


where  Fj  ■  (I  -  k8k8)  •  (k  A  nj),  g  »  k£  -  kg,  k  -  k^ A  e  and  fij  *  |z  “  dxdy/dsj  . 
Here  (x,y)  is  the  mean  level  plane  of  the  surface,  z  lying  orthogonal  to  (x,y).  This 
coordinate  system  is  centred  on  the  jth  element.  Further  k^  and  k8  are  the  incident 
and  scattered  wave  number  vectors,  nj  is  the  unit  normal  vector  of  the  jth  element,  e 
is  the  unit  vector  of  the  incident  field  (ez  its  z  component).  The  integral  is 

evaluated  over  that  portion  of  the  level  plane  corresponding  to  the  projection  of  the 
elemental  surface  onto  this  plane.  To  evaluate  this  integral  the  method  of  stationary 
phase11  can  be  used  if  k  is  large. 

*  a  is  the  backscatter  cross  section  for  a  unit  area  of  scattering  surface,  o  is 

defined  by  the  relationship  a  -  4irr2(E8>  /<E?>  where  Eg  is  the  scattered  electric 
field  intensity  derived  for  a  unit  area  of  the  scattering  surface  at  distance  r  from 
the  surface  and  E^  is  the  incident  electric  field  intensity  seen  by  this  unit  area. 
a  is  also  known  as  the  backscatter  coefficient. 
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Using  this  method  of  integration  the  solution  of  Ejg  is: 


E.  -  e.(e^kr/r)F./q  •  (n.  •  e  )  exp  ik(q  •  c.)|a.C.  -  B? 
js  3  “J r  -J  ~z  -  *]  I  3  3  3 


2  2 

where  e.  »  +1  at  a  relative  maximum,  and  A.  -  3  k./3x  ,* 
J  3  2.  3 

-1  at  a  relative  minimum,  B.  -  3  k./3x3y. 


-1  at  a  relative  minimum, 
+i  at  a  saddle  point. 


J  2  J  j 
C.  -  3  k./3y  . 
3  3 


To  determine  the  total  scattering  field  Eg  ,  the  contributions  from  each  scattering 
element  are  summed.  Consequently: 


(e*)  =  |F|2/q2r2(cos  ez)4 


E  iDiDi 


k!  exp(ikq  •  (r.  -  rk)) 


where  cos  0  =  q  •  e  as  F./n.*  e  =  F/q  •  e  and  D.  ■  A.C.  -  B. 

z  —  — z  J  —  J  — z  —  -Z  J  J  J  J 


Since 


o  =  4nr 


a  -  4n|F] 2yq2(cos  6 ^)4 


\  E  lDJD^rl  CXP(i-q  '  (^j  "  Sk)} 
W  i 


Assuming  the  phase  to  be  random  between  each  of  the  scatterers,  ie  each  scatters 
independently  of  the  others.  This  is  reasonable  as  the  phase  dispersion  introduced  by  a 
rough  scattering  surface  is  large  as  the  characteristic  scale  size  in  most  cases, 

and  over  a  large  number  of  scattering  elements  the  cross  terms  will  sum  to  zero.  Thus: 


a  *  4tt|f|  2^q2  cos462 


I  |DjV* 


Expressing  this  in  terms  of  the  two  principal  radii  of  curvature  at  each  point,  the  radar 
cross  section  becomes 


*(  E|(RaVi(R.V. 


*  kj  is  the  elevation  of  point  j  above  the  mean  level  plane. 


12 


where  R  /R.  -  1/K,  and  D.  cos  8  -  K.  |(3r/3x)  a  (3r/3y) |  ,  K.  being  the  Gaussian 

ah  J  z  J  j 

curvature  at  point  j  . 

Taking  the  ensemble  average: 


a  *  irRaRb<N> 

where  N  is  the  average  number  of  specular  points  per  unit  area.  Finally  introducing 
F(n)  the  Fresnel  reflection  coefficient,  then 

a  »  ^|Fj(n)  |2^iriTR^<N>  , 

n  being  the  angle  of  incidence  to  the  surface  for  each  point  j  . 

To  be  of  practical  use  and  (N>  must  be  re-expressed  in  terms  of  more 

easily  observable  quantities.  Work  by  Barrick1^  has  shown  that  these  factors  can  be 

4 

replaced  by  the  expression  sec  yp(Ax,Ay)  to  a  good  order  of  approximation.  Here  y  is 
the  angle  between  the  normal  to  the  level  plane  of  the  surface  and  the  local  surface 
normal  at  the  scattering  point.  For  backscatter  y  equals  the  angle  of  incidence.  The 
expression  p(Ax,Ay)  is  the  joint  slope  probability  distribution,  where  Ax  and  Ay  are 
the  surface  slopes  measured  in  orthogonal  directions  within  the  surface  at  some  point. 

The  expression  for  a  now  becomes: 

0  -  Tr^|F.(n)|  2>  sec\p(Ax,Ay) 

This  expression  can  be  used  to  determine  the  mean  waveform.  In  using  this  expression 
within  the  illuminated  area  the  wavefield  will  be  assumed  to  be  homogeneous.  As  o  is  a 
statistical  parameter  owing  to  its  dependence  on  p(Ax,Ay)  the  resulting  waveform  has  an 
averaged  form.  Finally  as  a  ■  Jc(z)dz  ,  z  being  orthogonal  to  the  plane  of  the  mean 
surface  then  cr(z)  will  take  the  form: 

^lFj(h)  |  it  sec\p(z;Ax,Ay)  . 

A  description  of  p(z;Ax,Ay)  must  now  be  obtained. 

3.3  Description  of  the  joint  height,  slope  probability  density  function 

In  past  work  p(z;Ax,Ay)  has  been  assumed  to  be  normal.  Insitu  measurements  of 

sea  surface  height  elevations  as  a  function  of  time  have  shown  the  wavheight  pdf 

13 

(probability  density  function)  to  have  a  skewed  Gaussian  form  .  Consequently  the  joint 
height  and  slope  pdf  (hereafter  known  as  wave  joint  pdf)  will  not  be  normal.  In  practice 
the  crests  and  troughs  of  a  rough  sea  surface  are  of  unequal  curvature.  The  troughs  are 
far  more  rounded  than  the  crests  which  often  exhibit  a  marked  degree  of  peakiness.  If 
the  waveheight  pdf  were  normal  the  crests  should  be  mirror  images  of  the  troughs.  To 
explain  this  observed  skewness  in  the  waveheight  pdf  non-linear  interactions  between  the 
wave  motions  interacting  to  produce  a  sea  surface  are  invoked. 

J  -i  - 
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Longuet-Higgins  has  derived  the  form  of  p(z;Ax,Ay)  for  a  long  crested  see  (here 
the  cross  section  of  a  sea  surface  is  taken  to  be  dependent  on  x  and  z  only, 
consequently  Ay  •  0)  and  finds  that: 

p(z;Ax)  -  (2lT02°^x)  j1  +  (*/6)  ((z/oz)3  -  3*/°z)  +  (X  j  z/2oz)  ((Ax/o^)  2  -  ljj 

X  exp^-  i((z /a2)  +  (Ax/aAx>2)j  • 


Here,  and  o^x  are  the  standard  deviations  of  the  surface  elevations  and  slopes  and 

X  and  X  the  skewness  coefficients  defined  by: 


M30/az  * 


X1  ‘  vl2/o*V02 


where  u  is  a  moment  of  the  distribution  p(z;Ax)  given  by  ( zmAx  )  . 
mn 

Jackson*5  has  used  the  results  of  Longuet-Higgins  with  the  assumption  of  a  fully 
developed  sea*  to  show  that  X  -  2X  .  For  radiation  to  be  received  back  at  the  altimeter 
0  <  0.5°  consequently  only  the  backscatter  component  is  of  interest  in  altimetry,  thus 
Ax  «  Ay  •  0  to  a  good  approximation.  So 

p(z;0)  «  (2 wo2 )  ^1  +  (X/6)  ^(z/oz)3  -  9z/ °z)^  exp(“  z2 /2o2)  • 

Therefore  the  final  form  for  o(z)  is: 

(|Fj(0)|2)  irp(z;0)  ,  as  sec2y  «  1  and  n  ■  0  . 

-  Oq (2it°2)  ^1  +  (X/6)  ^(z/oz)3  -  9z/ e*p(~  z2/2o2) 

This  above  pdf  has  its  mean  located  at  z  -  -Xaz5  *  while  the  mean  for  a  normally  dis¬ 
tributed  pdf  is  at  0m.  This  difference  in  the  position  of  the  mean  is  known  as  the 
wave-bias. 

3.4  Derivation  of  the  flat  sea  impulse  response 

Fig  6  illustrates  the  geometry  associated  with  the  reflection  of  a  pulse  of 
infinitesimal  duration  (strictly  a  mathematical  abstraction  produced  for  theoretical  con¬ 
venience)  by  a  flat  sea  surface.  0A  is  the  boresight  axis  of  the  altimeter,  ie  the  axis 
along  which  the  gain  of  the  altimeter's  antenna  is  a  maximum,  OB  the  axis  along  which 
power  contributed  by  element  B  arrives  and  Ef  is  the  attitude  error  angle.  Taking 
the  power  output  of  the  altimeter  transmitter  as  P  watts,  the  gain  of  the  antenna  as  G 
then  the  effective  radiated  power  is  PC  watts.  Finally  the  pulse  is  transmitted  at  time 

*  A  fully  developed  or  aroused  sea  is  one  whose  energy  spectrum  has  attained  equilibrium. 
Physically  this  means  that  a  balance  has  been  reached  between  the  energy  passed  into 
the  waves  by  the  wind  to  that  lost  when  the  waves  break. 
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zero.  In  deriving  this  response  the  method  used  here  is  a  simplification  of  that  used  by 

_  16 
Brown  . 

2 

At  B  the  power  received  per  unit  area  will  be  PG/4wR.  ,  the  power  reflected  back 

2  1 

from  an  area  dA  at  B  being  PGadA/4rR.  .  Thus  the  power  per  unit  area  received  back 

2  2 

at  the  altimeter  will  be  PGodA/4irR,  x  l/4irR,  .  For  an  aerial  of  gain  G  ,  working  at  a 

1  2  1 

wavelength  X  ,  the  effective  area  is  GX  /4ir  .  Consequently  the  power  received  through 
the  aerial  is  PG/4irR2  *  GX2/4it  *  cdA/4irR2  . 

For  such  an  infinitesimal  pulse  being  reflected  from  a  flat  sea  surface,  the  region 
of  the  sea  illuminated  is  initially  just  a  point  which  then  becomes  a  ring  of  infinitesi¬ 
mal  width  and  increasing  diameter.  The  area  of  such  a  ring  is  pd$  .  If  the  distance  of 
the  points  of  a  ring  from  the  altimeter  is  Rj  (Fig  6),  then  power  is  received  from  the 
ring  2Rj/c  seconds  after  the  pulse  is  transmitted.  Consequently  the  power  P(t) 
received  from  the  ring  at  time  t  is: 

P(t)  -  J<5(t  -  2R1/c)(pG2aX2/R^(4^)3)pd4. 


where  6(t  -  2Rj/c)  sets  the  time  of  reception,  since  S(t  -  2Rj/c)  is  unity  if 
t  =  2Rj/c  and  zero  otherwise.  From  the  geometry, 

R4  -  H4(l  +  (p/H)2)  . 

The  gain  of  the  aerial  is  given  by 

G(8)  -  Gq  exp(-  2  sin  6/y) 


where  6  is  the  angle  between  the  boresight  axis  and  the  OB  axis  in  Fig  6.  For  a  flat 
sea  and  small  attitude  error  (<0 . 5)  a  is  virtually  constant.  Setting  a  ••  to  denote 
this  and  substituting  K  ■  p/H  then 


P(t) 


2  2 
AG 

<4ir)3H 


SL  [ 

J 


<5(t  -  2Rj/c)  exp{-  4  sin2e/rj 

(1  +  K2)* 


pd$ 


Before  this  integral  can  be  evaluated  the  functional  dependence  between  B,  p  and  4 
must  be  established.  From  Fig  6: 


As 


also,  then 


AB2  -  p2  +  H2  t«n2ER  -  2pH  tan  ER  cos  *  . 


AB2  -  R2  +  R2  -  2R,R2  cos  8 


cos 


8  -  ^  +  R2  -  p2  -  H2  tan^jj  ♦  2pH  tan  ER  cos  ♦jj^lRjRj  . 
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But 


and  so 


R2  -  H  sec  Er  and  R2  -  H2(i  +  K2) 
:os  6  -  ^cos  ER  +  K  sin  ER  cos  1  +  K2)  ^ 


Hence  the  resulting  integral  is 


P(t) 


X2G2aQ  ^  6(t  -  2Rj/c)P  exp{w| 


(4s) 


,a0  pu  -  ZKj/cit'  t 

j  o  +  kV 


pd<J> 


where 


-  -  < 

y 


,  J 


cos  ER  +  K  sin  ER  cos  <|ij 


1  +  K 


This  expression  is  known  as  the  impulse  response  integral.  As  shown  in  Appendix  A: 


P(t) 


— ^  exp.  -  ±  sin2ERj  expj-  ^  cos2ERJ  •  (1  +  a2) 


2  2 

where  a  »  4ct  sin  2Er/y  H  . 

This  description  is  valid  for  t  >  0  .  For  t  <  0  ,  P(t)  »  0  . 


3.5  Final  convolution  procedure 

Having  obtained  the  quantitative  description  of  the  impulse  response  three  further 
convolutions  must  be  performed  before  a  quantitative  description  of  a  mean  waveform 
P(t)  is  obtained.  Schematically  P(t)  is  given  by 


P(x)  -  Ir(T)  *  Ps(t)  *  W(t)  *  1 1 ( t)  , 

where  Ir(T)  is  £he  flat  sea  impulse  response,  Pg(x)  is  the  pulse  shape  distribution, 
W(t)  is  the  waveheight  pdf  and  It(x)  is  the  instrument  transfer  function.  The 
symbol  *  denotes  a  convolution  process.  The  order  in  which  the  above  convolutions  are 
carried  out  is  not  important  on  physical  grounds,  however  mathematical  simplicity  does 
dictate  an  order.  This  is  that  firstly  the  convolution  of  Pg  with  W  is  performed, 
this  produces  a  function  R  .  Then  R  is  convolved  with  I  producing  a  function  S 
and  lastly  S  is  convolved  with  giving  P  .  The  formation  of  R  is  given  in 

Appendix  B,  that  of  S  and  P  in  Appendix  C. 

Two  points  must  be  noted.  In  approaching  a  roughened  sea  surface,  surface  eleva¬ 
tions  represented  by  the  tail  on  the  positive  side  of  the  time  axis  of  a  waveheight  pdf 
will  be  encountered  first  by  an  impinging  altimeter  pulse.  As  the  convolution  process 
causes  a  function  to  be  reversed  before  integration  proceeds,  for  P  to  portray  the 
effect  of  an  asynmetric  pdf  correctly  the  waveheight  pdf  must  be  reflected,  it  reversed 
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about  its  probability  axis  before  the  convolution  takes  place.  In  Appendix  C  this  is  done 

by  multiplying  the  parameters  F  and  G  by  minus  one.  Secondly  the  function  describing 

P  has  been  described  by  a  normal  distribution.  Consequently  it  is  symsetrical. 

8 

4  ACCURACY  OF  THE  SEASAT  1  ALTIMETER  HEIGHT  MEASUREMENT 

4.1  General 

The  Seasat  I  processed  altimeter  height  data,  ie  the  initial  height  measurement 
corrected  for  instrumental  and  ionospheric  and  tropospheric  propagation  delays  is  not  a 
measurement  of  the  true  satellite  altitude,  ie  distance  from  the  spacecraft  CG  to  mean  sea 
level.  Two  further  corrections,  the  tracking  and  wave-bias  corrections  have  to  be  applied 
to  the  processed  height  values  to  derive  true  height  values.  Consequently  the  accuracy  of 
the  true  altitude  measurement  depends  on  the  accuracy  with  which  each  of  these  corrections 
is  known,  this  accuracy  being  expressed  in  terms  of  a  standard  deviation  in  the  two 
corrections. 

The  ideal  tracking  situation  would  be  for  the  ATU  to  automatically  track  to  that 
point  on  the  leading  edge  of  a  waveform  hereafter  known  as  the  mean  sea  point  whose 
associated  power  level  P  corresponds  to  the  contribution  of  those  specular  points 
(.facets)  being  at  mean  sea  level.  Unfortunately  the  time  ordinate  of  this  point  varies 
with  the  prevailing  sea  state  conditions  (ow  and  X)*  and  the  attitude  error  angle  ER  . 

So  without  knowledge  of  these  three  variables  it  cannot  be  tracked.  In  setting  up  the 
altimeter  on  Seasat  I  the  variables  X  and  ER  were  taken  to  be  zero  and  by  using  AGC 
to  keep  the  mean  integrated  power  per  100  consecutively  received  pulses  constant  the 
point  on  the  leading  edge  of  a  mean  waveform  whose  normalised  power  level  was  0.443  was 
tracked.  The  processed  height  measurements  were  corrected  during  ground  based  processing 
for  the  effect  of  ER  i  0°  but  X  was  always  assumed  to  be  zero  as  the  waveheight  pdf 
was  always  assumed  to  be  normal. 

The  work  presented  here  although  similar  to  that  carried  out  by  Hayne * ^  was 
developed  independently.  The  derivation  of  the  quantitative  description  of  a  mean  wave¬ 
form  is  presented  here  in  greater  detail  than  in  the  published  work  of  Hayne.  Further  a 
description  of  the  discriminator  curve  is  also  given. 

If  the  ATU  is  arranged  to  track  the  mean  sea  point  for  a  particular  value  of 
where  E_  -  0°  and  X  ■  0°  then  for  other  values  of  o  or  either  E_  J1  0°  or  X  +  0 

K  W  R 

for  the  same  value  a  different  point  will  be  tracked.  In  each  case  the  time 
separation  between  the  point  tracked  and  the  mean  sea  point  must  be  calculated.  This 
distance  equals  the  sum  of  the  tracking  and  wave-bias  corrections. 

4.2  Tracking  and  wave-bias  corrections 

As  mentioned  in  section  2.4  the  tracking  law  employed  by  the  ATU  is  of  the  form 
I(t)  ■  M(t)  -  oL(t  ♦  W)  I(t)  is  the  output  of  the  tracking  gate  for  time  offset  t  , 
t  being  the  displacement  of  the  centre  of  the  middle  gate  M  from  the  origin  of  the 
time  axis  of  the  waveform.  I.  is  the  late  gate  output,  W  the  width  of  the  late  gate 
and  a  a  parameter  which  in  conjunction  with  W  dictates  which  point  on  the  leading  edge 

*  X  here  denotes  the  waveheight  pdf  skewness  variable. 
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is  tracked.  The  distribution  of  I  with  t  is  known  as  the  discriminator  curve  and  is 
obtained  by  convolving  the  waveform  with  the  middle  and  late  gate  functions  and  then 
applying  the  expression  for  I(t)  given  above.  The  deviation  of  the  expression  of  a 
discriminator  curve  is  discussed  in  Appendix  D.  Both  the  waveform  and  discriminator 
curve  have  coincident  time  origins.  The  tracking-bias  is  given  by  the  value  of  t  when 
I(t)  »  0  ,  with  the  point  on  the  waveform  corresponding  to  I(t)  -  0  being  tracked. 

This  bias  value  is  dependent  on  the  values  of  X,  a  ,  ER,  W  and  a  . 

The  wave-bias  correction  arises  due  to  the  displacement  of  the  mean  of  a  skewed 

normal  distribution  compared  to  the  mean  of  a  normal  distribution.  If  the  waveheight  pdf 

is  represented  by  a  normal  distribution  then  the  point  on  a  waveform  whose  time  ordinate 

is  at  the  origin  of  the  time  axis  is  the  mean  sea  point.  For  a  waveheight  pdf  having  a 

skewed  form  the  time  ordinate  of  the  mean  sea  point  will  be  2Xo  /c  seconds.  This  is 

w 

the  magnitude  of  the  wave-bias  correction. 

4.3  Simulating  the  determination  of  and  X 

To  apply  the  tracking  and  wave-bias  corrections  X,  and  ER  must  be  known. 

E  is  calculated  from  knowledge  of  a  spacecraft's  pitch,  roll  and  yaw  angles  as  measured 
by  the  spacecraft's  attitude  stabilisation  system.  The  values  of  X  and  and  their 

variances  are  obtained  by  a  least  squares  fit  in  these  parameters  to  mean  waveform  data. 
To  answer  the  question  as  to  the  accuracy  of  a  true  height  value  obtained  from  processing 
Seasat  1  altimeter  data  this  fitting  process  was  simulated  to  provide  values  of  the 
standard  deviation  of  the  tracking  and  wave-bias  corrections. 


I 
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Examination  of  the  form  of  a  mean  waveform's  quantitative  description  given  in 
Appendix  C  shows  that  is  is  non-linear  in  the  variables  X  and  a w  .  Consequently  this 
necessitates  use  of  the  differential  least  squares  procedure  so  that  the  normal  equations 
developed  in  this  method  are  linearised,  ie  only  linear  terms  in  X  and  o w  are  produced. 
However  use  of  the  differential  least  squares  method  does  not  automatically  guarantee 
that  sensible  or  realistic  values  of  the  fitting  variables  will  be  determined  as  the 
procedure  is  iterative  and  can  thus  diverge.  Although  only  a and  X  are  of  interest, 
to  fit  the  theoretical  waveform  to  modelled  waveforms  a  third  variable  T  (time  displace¬ 
ment)  must  be  introduced  as  the  origin  of  the  modelled  waveform's  time  axis  is  not  known. 

If  uncertainty  occurs  in  the  integrated  power  level  a  fourth  variable,  power  gain  must  be 
fitted  for.  Adopting  this  method  a  matrix  equation  of  the  form  DX  »  R  results,  where 


YY 

Y"1  3f  3f 
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j 


£  is  the  function  describing  the  theoretical  waveform  and  3f/3ow,  3f/3T  and  3f/3X 

the  derivatives  at  points  i,  i  -  I,  60  (i  ■  I  corresponds  to  gate  -30,  i  •  60  to  gate 

+30  of  the  DFB  of  the  Seasat  1  altimeter  system) .  Ac  ,  AX  and  AT  are  increments  in 

w 

the  guess  values  of  o^,  X  and  T  to  be  solved  for,  and  r  are  the  residuals. 

To  determine  the  standard  deviation  and  correlation  between  the  fitted  variables 
a  covariance  matrix  was  introduced.  As 


D-1  DX  -  D-1  R  . 

Consequently  X  «  (D  1  D)  1  D  1 R  ,  the  covariance  matrix  C  being  (D  1  D)  * .  Having 

found  the  values  of  a  ,  X  and  T  the  variance  of  each  of  these  terms  is  found  by  first 

w 

calculating  the  variance  of  the  residuals  ,  ie  using  the  above  o^,  X  and  T  values 
in  the  theoretical  description  and  determining  the  residuals  between  the  predictions  of 
this  description  and  the  modelled  data. 

Writing  C  as 


II 

*12 

‘21 

*22 

l31 

*32 

,  of 

AX 

is 

A°w 

and 

AX 

and 

AT 

is 

then  the  variance  of  Ao^  is  vr*j|  >  of  AX  is  vr®22  *  *nd  AT  **  Vr*33  * 

The  correlation  coefficient  between  Acw  and  AX  is  a|  la22  *  between  A<JW  and 
AT  is  *13/^*11*33  *nd  between  AX  and  AT  is  *23/^*22*33  * 

To  start  the  fitting  procedure  initial  guess  values  had  to  be  provided.  Initially 

no  noise  was  added  to  the  modelled  waveform  but  even  in  this  case  the  method  diverged  - 

the  normal  equations  were  virtually  multiples  of  one  another  as  the  determinant  derived 
•  1 

in  producing  the  inverse  of  D  D  became  singular  after  a  few  iterations.  Only  in  the 
case  of  0v,  X  and  T  being  very  close  to  the  values  used  in  producing  the  noise-free 
waveform  did  the  method  converge  and  produce  the  correct  values  of  o^,  X  and  T  .  To 
overcome  this  singularity  difficulty  a  weighting  matrix  was  introduced. 
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As  D  X  *  R ,  let  W  be  a  weighting  matrix  then: 

^  -"w-w  »m>  'V 

HDX  -  UR 

and 

■V  «V—  I  *v<v  «  1  1  'V'V 

X  -  (D  WD)  D  WR 

•s,  |  <v<v  _  |  ^ 

the  covariance  matrix  C  now  being  (D  WD)  •  W  is  a  diagonal  matrix  in  this  case 
having  dimensions  3  *  3,  After  some  initial  experimentation  it  was  found  that  provided 
Aow  and  AX  were  highly  weighted  (implying  little  change  in  their  value  over  successive 
iterations)  compared  to  T  convergence  was  possible.  In  practice  the  weighting  factors 
applied  to  Ac^  and  AX  were  equal  having  the  value  1 0^ ,  while  the  third  weighting 
factor  had  a  value  of  unity.  Varying  these  values  did  not  affect  the  value  of  the 
variances  in  Ao^,  AX  and  AT  significantly  provided  the  ratio  of  the  weighting  factors 
of  Aow  and  AX  to  AT  were  kept  >10^.  As  this  ratio  was  reduced  singularities  were 
more  likely  to  occur  with  the  iterative  process  often  diverging  beforehand.  Random  noise 
was  then  added  to  a  modelled  waveform  -  as  in  practice  waveforms  are  contaminated  with 
noise.  The  independence  of  the  variance  of  Ao^,  AX  and  AT  with  the  weighting  factor 
values  was  then  verified.  Further  the  values  of  Ao^,  AX  and  AT  were  all  highly 
correlated  owing  to  the  use  of  such  large  weighting  factor  values  applied  to  two  of  the 
variables.  The  values  of  the  variances  however  were  unrealistic  as  the  variance  on  each 
of  the  variables  was  much  smaller  than  the  difference  between  the  value  of  the  variable 
used  in  deriving  the  model  waveform  and  the  value  obtained  from  the  fitting  process. 

To  derive  more  realistic  values  of  the  variances  the  following  method  was  employed. 

For  specified  values  of  o^,  X  and  T  a  modelled  waveform  was  produced.  Random 

noise  was  added  to  this  waveform  to  attain  a  SNR  of  either  30,  20  or  10.  Initial  guess 

values  were  derived  by  adding  random  noise  to  the  modelled  waveform's  values  of  o^,  X 

and  T  .  For  this  guess  value  was  within  I  OX  of  the  model  value,  for  X  within  50X 

and  T  within  I00X.  With  these  guess  values,  convergence  was  obtained  after  10  to  15 

iterations.  This  process  was  then  repeated  a  further  99  times,  using  the  same  o^,  X 

and  T  values  but  using  different  initial  guess  values  and  adding  new  noise  to  the 

modelled  waveform  producing  the  same  SNR  as  before.  The  result  of  this  process  was 

100  o^,  X  and  T  values.  Having  such  a  large  data  base  enables  reliable  values  of 

the  variance  in  o  ,  X  and  T  to  be  calculated, 
w’ 

Having  done  this  for  one  set  of  and  X  values  this  procedure  was  repeated 
for  further  values  of  o^  ,  Owing  to  the  high  correlation  between  Aow  and  AX  hence 
o  and  X  ,  the  tracking  and  wave-bias  corrections  will  be  highly  correlated.  Conse¬ 
quently  the  mean  and  variance  of  the  wave-bias  and  tracking  correction  must  be  obtained 
by  calculating  the  two  corrections  for  each  of  the  and  X  values  and  then  finding 
the  mean  and  variance  of  the  biases  rather  than  finding  the  mean  and  variance  of  both 
ow  and  X  from  the  100  values  and  then  determining  the  mean  and  variance  of  the  two 
biases.  Similarly  in  calculating  the  mean  of  the  wave-bias  plus  tracking  correction, 
ie  the  offset  correction,  this  correction  must  be  determined  for  each  ow  and  X  value 
and  then  the  mean  and  variance  calculated. 
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One  further  possibility  remains,  does  the  magnitude  of  the  standard  deviation  of  the 

offset  correction  vary  significantly  with  X  for  a  constant  value  of  a  and  E  ? 

■  **  © 

To  answer  this  question  this  standard  deviation  was  calculated  for  X  -  0  and  ER  ■  0  . 

The  results  obtained  for  the  standard  deviation  for  X  -  0  and  0.4  were  similar 

(differing  by  only  a  few  centimetres)  to  the  results  obtained  for  X  -  0.2  .  To  provide 

estimates  of  the  error  bars  on  the  points  in  a  plot  of  rms  waveheight  versus  standard 

deviation  of  the  offset  correction  the  results  for  SNR  -  30  and  20  have  been  used  to 

predict  results  at  SNR  -  10  -  this  being  the  mean  SNR  for  Seasat  1  mean  waveforms.  For 

SNR  »  30  the  values  of  the  standard  deviation  of  the  offset  correction  were  multiplied  by 

/T  ,  for  SNR  •  20  by  /J.  The  results  for  X  *  0,  0.2  and  0.4  with  SNR  -  10  were  then  used 

to  derive  mean  values  of  the  standard  deviation  of  the  offset  correction  and  its  standard 

deviation  for  various  o  values. 

w 

4.4  Tracking  error  contribution 


To  determine  the  accuracy  of  the  true  height  measurement  one  further  error  source 

must  be  considered  -  tracking  error.  This  arises  owing  to  the  presence  of  noise  on  the 

received  waveform.  This  noise  causes  an  error  in  the  positioning  of  the  tracking  gate, 

1 8 

Le  value  of  t  for  the  condition  I(r)  >0  .  From  the  work  of  Brooks  and  Dooley  the 
standard  deviation  of  this  height  error  o^  introduced  by  random  noise  is: 


T  At 
e 

12 


metres. 


where  a  is  the  power  SNR,  Tfi  and  T.  are  the  widths  of  the  middle  and  late  gates 
expressed  in  metres.  At’  is  the  height  resolution  cell  size,  the  rms  waveheight, 

a  is  a  constant  of  value  0.3327,  T^  is  the  point  tracked  to,  and  N^  is  the  number 
of  pulses  averaged. 


In  this  work 
-8 


T  has  been  taken  to  be  0.5.  At  •  c  where  t  *  1  ns  ,  and 
P  t 

3  x  10”  metres.  T  •  T.  -  a.c  where  a.  are  the  number  of  DFB  gates  combined  to 
e  l  i  t  l 

form  the  middle  and  late  tracking  gates  and  N  -  1000 

P 


10  °w  / 0.5  SNR  +  \\  aj(ct)' 

(0.3327)2  \  535  /  \ai/  r5~" 


With  these  values: 

I 

metres. 


The  values  of  for  various  values  of 

these  values  is  discussed  in  Appendix  D. 

4.5  Determination  of  the  total  error 


ow  are  given  in  Table  2.  The  reason  for 


Finally  to  arrive  at  the  true  height  accuracy  the  variances  of  the  offset  correc¬ 
tion  and  tracking  error  were  added,  the  square  root  of  the  sum  being  taken  to  give  a 
figure  for  the  above  accuracy.  The  assumption  made  here  was  that  the  offset  correction 
and  tracking  error  were  statistically  independent.  Table  3  indicates  the  figures  of 
accuracy  expected  for  a  variety  of  o^  values  with  SNR  -  10.  These  results  are  plotted 
in  Fig  7. 
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This  figure  shows  that  Che  true  altitudes  in  the  case  of  the  Seasat  1  altimeter  can 
be  calculated  with  an  accuracy  of  better  than  20  cm  provided  ow  <  5  m.  Consequently  it 
is  worth  applying  the  total  correction  as  the  magnitude  of  this  is  greater  than  the  magni¬ 
tude  of  its  standard  deviation.  However  it  must  be  remembered  that  Fig  7  has  been 
obtained  assuming  that:  (1)  the  transmitted  pulse  shape  is  normally  distributed,  (2)  that 
the  waveheight  pdf  used  in  the  waveform  derivation  is  valid  for  all  sea  state  conditions, 
and  (3)  that  the  theoretical  treatment  of  the  backscattering  problem  used  here  is  adequate. 
If  the  pulse  shape  of  (1)  were  not  normal,  no  effect  on  Fig  7  is  to  be  expected.  To 
verify  points  (2)  and  (3)  further  experimental  work  is  required. 

5  CONCLUSION 

This  Report  has  presented  the  development  of  a  more  realistic  description  of  a  mean 
waveform  than  those  developed  previously.  Having  done  this  the  question  -  to  what 
accuracy  can  the  true  altitude  be  calculated  from  Seasat  I  altimeter  height  measurements? 
has  been  answered.  For  the  SNR  applicable  to  this  altimeter's  mean  waveforms  it  is 
found  that  the  accuracy  is  better  than  10  cm  for  o  »  1  m,  while  even  for  very  rough  seas 
(a  %  5  m)  the  accuracy  is  better  than  20  cm. 
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Appendix  A 

EVALUATION  OF  THE  IMPULSE  RESPONSE  INTEGRAL 

In  section  3.4  a  quantitative  description  of  the  impulse  response  is  desired.  To 
obtain  this  an  integral  of  the  form 


P(t) 


iVo  f  Pa, 

(4ir)  V  J  (  1  +  Kj  L 

P 


has  to  be  evaluated.  This  Appendix  is  concerned  with  this  evaluation.  In  the  above 
expression: 

.2, 
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1  - 


(cos  ER  +  K  sin  ER  cos  d>) 
I  +  K2 


1  +  K2. 


4K  sin  2E_  cos  4K2  sin2E  cos2it 

+  - L_ -  +  - £— - 


Yd  ♦  K  ) 


Yd  +  K  ) 


where  K  =  p/H  . 

A  point  at  distance  R  from  the  altimeter  will  return  a  contribution  at  time 
t  -  2R/c  seconds  after  pulse  transmission. 

As  R  -  (H2  +  p2) *  and  K  -  p/H  then 


Translating  the  time  origin  such  that  time  measurements  are  now  made  with  respect  to  the 
time  at  which  the  mean  sea  level  contribution  is  received,  ie  t  ■  t  -  2H/c  ,  H  being 
the  distance  between  the  spacecraft  centre  of  gravity  and  the  mean  sea  level,  then 
K  -  (c2t2  +  4Htc)^/2H  .  As  H  »  8  x  105  m  (Seasat  1  satellite)  and  |t|  <  100  ns 
(width  of  Seasat  1  DFB  dictates  the  time  range  in  which  a  signal  can  be  sailed)  then 

i  2  4 

K  »  (cx/H)5  to  a  good  order  of  approximation.  Further  K/(l  ♦  K  )  *  (cx/H)’  and 
K2/ ( 1  +  K2)  -  cx/H  .  Consequently 


P(x) 


2  2 

.  ,  xVo 

where  A  *  ■  v  ;  exp 


(4ir)  V 


'  /  «( x  -  |  (R  -  H))  fexp  sin  2ER  cos 

$ 

*  expj^  sin^  cos  2f|  pdf  , 

{"  7  8in\}  expf  IT  co>2er}  ’ 


.2,2 


(I  +  K  )  >1  to  a  good  approximation. 
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This  integral  can  be  evaluated  simply  by  expressing  each  of  the  exponential  terms 
as  a  power  series  and  integrating  the  terms  independently.  Taking  the  Seasat  I  satellite's 
orbital  and  altimeter  parameters  to  be  the  nominal  one,  the  maxi  tm  an  value  of 
(4/y) (ct/H)  sin2E„  is  ~2  *  10~6  where  4/y  -  7.112  x  103,  E»^  ■  i°.  c  -  3  x  108, 
t  -  10-/  seconds  and  H  ■  8  *  10  .  Consequently  the  contribution  of  exp((4/y) (ct/H)  x 

1DAX 

sin2E„  cos  p)  will  be  virtually  constant,  ie  independent  of  t  and  have  a  value  of  unity 
R 

to  a  very  good  order  of  approximation.  Thus  the  integral  reduces  to: 


P(t)  -  A 


Js  (T  “  "c  (R  ”  H^)  eXp|  yfir)  8in  2ER  008  pd<P  * 


A  /ct  \a 

Let  Oj  -  am  2ER  ,  then: 


As 


2tt 

f  n. 
f  cos  p 


2  „3 

al  2  al  3 

exp(aj  cos  p)  m  1  +  Oj  cos  <p  +  -jj-  cos  4>  +  yr  cos  I1  +  •  •• 


dip  is  zero  if  n  is  odd,  and 


2tt 

2ir 

2  IT 

/  cos^$  d$  -  it  , 

f  4  3 

/  cos  p  dp  -  J  V  , 

J 

f  cos 8p  dp  ■  ■j  i 

J 

0 

0 

J 

0 

and 


then 


2ir 

/ 


3  35 

cos  $  dp  -  64  ^  * 


2ir  2  4  6  8 

f  al  al  3  al  5  °1  35 

J  exp (a,  cos  P) dp  -  + 


Consequently 


P(t)  ■  2ft A  ^1  ♦  a2  +  yg  *  yyg  +  .  . 


Act  2 

where  a  -  —5—  sin  2E  . 

A  R 

As  the  maximum  value  of  a  •  0.38  for  |t|  <  100  ns  and  the  preceding  Seasat  1 

values  then  a2  m  0.1451,  a8/ 16  •  0.0013  and  a8/576  ■  5  x  10  8  .  Thus  a2  >  a*  ►  a8 
2 

so  P(t)  *  2itA(  1  ♦  a  )  to  a  good  order  of  approximation. 


So  the  evaluation  of  P(t)  gives 


Appendix  A 


2lrX  G0C0 

- —  exp) 

(4  *>V  * 


-  sinV 

Y  R 


This  description  is  valid  for  t  >  0  .  For 


exp  -  cos ^  j(l  +  a2) 
t  <  0,  P(t)  -  0  . 
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EVALUATION  OF  THE  CONVOLUTION  OF  THE  WAVE  JOINT  HEIGHT  AND 
SLOPE~~PROBABILITY  DISTRIBUTION  AND  THE  TRANSMITTED  PULSE  DISTRIBUTION 

In  the  derivation  of  a  received  pulse's  mean  envelope  two  convolution  integrals  have 
to  be  evaluated  (see  section  3.5),  this  Appendix  is  concerned  with  the  first  of  these. 

The  two  functions  being  convolved  are  the  waveheight  pdf  (probability  density  function) 
and  a  function  describing  a  transmitted  pulse's  shape. 

From  section  3.3  the  waveheight  pdf  expressed  in  terms  of  time  is: 


(2™w)  exprS 


while  the  mean  pulse  shape  is  described  by  the  function: 


f  } 

((ct)3/8a2)-  (9ct)/2cw 


,  as  z  =  ct/2, 


(2%)  exP  -  77  • 


If  A  and  B  are  two  functions  that  are  being  convolved  then  from  Fourier 


transform  theory: 


FT (A  *  B)  =  FT (A)  *  FT(B) 


A  *  B  -  FT  (FT (A)  x  FT(B) ) 


Here  *  denotes  the  convolution  process,  FT  the  operation  of  a  Fourier  transform, 
and  FT  *  the  operation  of  an  inverse  Fourier  transform.  From  contour  integration  it  is 
found  that 

FT(exp(-St2))  =  (■£)  exp^-  , 


consequently 


2  2 

In  solving  this  last  integral  let  Q  -  c  /8ow  .  Rewriting  this  integral  in  its 
component  parts  one  gets: 
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(2lT°w) 


-i 


oo 

/ 


2  ^}c  /  3  2 

exp(-  Qt  )  exp(ju>t)dt  +  - s-  f  t  exp(-  Qt  )  exp(jwt)dt 

48cr  J 


^lc  f  2 

- |  t  exp(-  Qt  )  exp( jut)dt 


As 


f  2 

-A  (  »2\ 

I  exp(-  Qt  )  exp(jut)dt  ■  j 

q)  eXp(“  4Q/ 

(B—  1 ) 


solutions  to  the  integrals  involving  t  and  t  terms  can  be  easily  derived  as  follows: 
Letting 


<X} 

! 


exp(-  Qt  )  exp(jwt)dt 


then: 


d^  C  “  Jtc  ’ 


d2c  2 

.  _  t  c 

dw 


and 


d3c 


dw 


3  -  “jt3c 


.  I  2 

As  c  is  also  equivalent  to  (ir/Q)J  exp(-  w  /4Q)  then: 


and 


Jtc 


t2c 


jt3c 


"  GqXq)  exp(~  4q)  ’ 

' ')  exp("  fe) 


i 


Consequently  the  solution  of  integral  (B— 1 )  is: 

j3X.cw 


I  Si) 


I  - 


8o 


Hence  the  product  of  the  two  Fourier  transforms  is: 

2 


(I)  exp(~  (zij +  *2) 


j3X j  cw  X,c 


/  3w  w3  \ 

V*’7/ 
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Consequently  the  inverse  Fourier  transform  is: 


_1_ 

2ir 


oc 

/ 


U  exp(- 


Vw2)0 


+  jXw  -  jYw 


exp(-  juit)dw 


(B-2) 


where  0  -  2/c,  V  -  (1/4Q)  +  (ap2/2),  X  -  3XiC/8owQ(c2/24o2Q  -  l)  and  'i  -  XjC3/384Q3o; 
This  integral  is  solved  in  the  same  manner  as  the  previous  one.  Setting 


then 


and 


C 


l  exp(-  Vw2) 


exp (-  j<»t)dt 


H  -  -  -  -  (w)(?)i  exp(-  Iv)  ’ 

§  -  -  “2c  -  (£X?)*(w  - ')  exp(-  h) 

±-|  =  +  jiA  =  (lv)(f)(B"^2)eXP("^v)  ' 


Consequently  the  solution  of  (B-2)  is 


u 

2(Vtt)* 


Thus  a  function  describing  the  result  of  the  convolution  of  the  wave  pdf  and  pulse 
shape  distributions  has  been  produced.  This  result  is  used  in  Appendix  C. 
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DEVELOPMENT  OF  THE  QUANTITATIVE  DESCRIPTION  OF  A  MEAN  WAVEFORM 


As  mentioned  in  section  3.5  the  quantitative  description  of  a  waveform  Pw  results 


from  the  following  convolution  procedure: 


P  =  I  *  W.  *  P_  *  I_ 
w  p  h  D  T 

where  I  is  the  impulse  response,  W,  the  joint  wave  height/slope  pdf,  V  the  distri- 
p  n  u 

bution  describing  the  pulse  shape,  1^,  the  instrument  transfer  function,  and  *  denotes 

a  convolution  process.  The  convolutions  can  take  place  in  any  order.  The  convolution  of 

W,  with  P  is  given  in  Appendix  B;  this  Appendix  is  concerned  with  the  remaining  two 
h  D 

convolutions,  zc  I  *  T  *  I„  where  T  =  W,  *  P_  . 

p  T  h  D 

Evaluation  of  T  *  I_ 


From  Appendix  B: 


T  = 


2(Vir) 


exp 


("  Iv)  (‘  +  2V  +  2v(4Tv2  “  2v)j 


and 


T]  <  t  c  t2 


T  <  V  T  >  V 


Consequently  the  integral  representation  of  the  convolution  is: 


|  2 

J  ( 1  +  F(t  -  t)  +  G(t  -  x)3j  exp(-  e(t  -  x)2)dx 


where  F  =  (X/2V)  -  (3Y/4V2),  G  =  Y/8V3  and  e  -  1/4V  . 

Substituting  T  =  e^(t  -  x)  then  -  dT/e^  «  dx,  T^  *  e^(t  -  Xj)  and 

T,  "  ~  V  • 

So  the  integral  can  now  be  expressed  as: 


1 

7 


r 

LT2 


exp(-  T  )dT  + 


7 1 


T  T 

T  exp(- T2)dT  *  -J  I  T3  exp(-T2)dT 


The  solution  of  the  first  integral  is  in  terms  of  the  error  function,  ie 


1 

7 


1 1 

j  exp(-  T2)dT  -  0.5^J  (erf  (Tj)  -  erfttp)  . 
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For  the  second  integral 


T  exp(“  T  )dT 


Xjexp(-T^)-exp(-T2)J  . 


To  solve  the  third,  consider 

J-  T2  exp(-  T2) 
dT 


then 


therefore 


T3  exp(-T2)dT 


/ 

T2 

f  '  T3  exp(-T2)dT 

P  * 


2T3  exp (-  T2)  +  2T  exp(-  T  ) 


2  2 

.  2  f  T  exp(-  T2)dT  -  i  exp(-  T  ) 


4  jO  +  exp(_  TD  ■  0 +  t5)  exp(~  T>)! 


Consequently  the  solution  of  T  *  1^,  is: 

0.5^|  jerf (T,)  -  erf(T2>J  +  *  4  (’  +  *2))  6XP("  *2) 

-(£♦7  (l  +  T?))  exp(" T?)  * 

The  desired  quantitative  waveform  description  can  now  be  obtained  by  evaluating: 

sf  (I  ♦  a,t)  «*p<-  ft)  jyJJ  (erf (Tj)  -  erf(T'))  ♦  ♦  4)  (eXP("  ^  "  ***('  ^  j 

0 

»(-  »J2)  -  Ti2  exp(‘  Ti2)jdt  • 


+  ^  t*2  exp( 
e 


where  B 


and 


2'»2°0°0  « 

exp| 

(4h)3H^  2(:iV)i 

4c  2 

W  coa  E*  ' 

aJ 

4c  .  2 

-  -a-  sin 

yTi 

eVwj  -  t) 

(w, 

•  S  -  T|) 

e^(w2  -  t) 

(v2 

"  *  ~  X2) 
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This  integral  splits  into  six  parts: 


Part  (1) 


Part  (2) 


Part  (3) 


00 

I  Jj  J  (erf(T|)  -  erf(Tp)  exp(-  £t)dt 
0 

‘(i  *  4)  /  (ext(-  t;2)  •  M-  Ti2)) 
4  /  (t22  «**(-  n2)  -  Ti2  M-  t;2)) 


exp(-  f t)dt 


exp(-  ft)dt 


Part  (4)  J-Jj  a l  j  c(erf(Tp  "  erf(T'))  exp(-  ft)dt 
0 


Part  (5) 


Part  (6) 


B  (ll  +  4)  ai  /  c(exp("  T22)  “  exp("  Ti2))  exp(“  £t)dt 

J  t(r’22  exp(-  T^2)  -  t;2  exp(-  TJ2)^  exp(-  £t)dt  . 


Evaluating  part  (1) 


Integral  is  of  the  form 


j  erf (T) 
0 


exp(-  ft)dt 


As  T 


e^(w  -  t)  then  dt  »  -dT/e 


i 


and 


00 

/ 


erf(T)  exp(-  ft)dt 


1 

7 


j. 


erf  (T)  expl 


(-  ((.  -  j)) 


dT 


Set  y  •  fw,  n  =  f/e^  . 

As 


^  erf (T)  exp(-  (y  -  nT)) 


f 


erf (T)n  exp(-  (y  -  nT))  + 


2 

77 


exp(- 


T2) 


exp(-  (y 


nT)) 
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— OO 

erf (T)  exp(-  (Y  -  nT))dT  -  ^  “  77  j  exp(-  T2)  exp(-  (y  -  nT))dT 


♦  erf(T)  exp(-  (y  ~  nT))| 


Evaluating  j  exp(-  T2)  exp(-  (Y  -  nT))dT  . 


Completing  the  square  this  integral  can  be  written  as 


2  —CD  —OO 

kp(t  -  y)  /  exp(~  (T  ~  7)  )dT  ■  exp(r " Y)  /  exp(~ s2)ds 


e2w-(n/2) 


-  4  «p(t  '  Y)('  +  erf(eiw  ‘  l)) 


where  S  -  T  -  ~  .  Consequently 


erf(T)  exp(-  ft)dt  «  — j-  erf(elw)  -  exp 

ne!  L 


(t  ~  Y)('  ♦  erf  (e*w  ■  7))] 


Hence  part  (1)  equals: 


-  erf 


Evaluating  pert  (2) 


(eJw,)  -  exp^  -  fv^l  *  erf^w,  - 
(X2)  ♦  exp^  -  fv2)^l  ♦  erf^e5w2  - 


This  integral  is  of  the  form 


exp(-  T  )  exp(-  f t)dt  . 


As  T  *  e^(w  -  t)  by  completing  the  square  this  integral  becomes: 
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exp(-  e(t  +  a) 


2 


)dt 


where  a  -  -  2w)  . 

Substituting  S  -  e^(t  +  a)  one  gets 

OO  2 

-y  exp^  -  Y)  f  exp(-  S2)dS  -  exp^  -  yjjjj  O  -  erf(e*a))  . 
e^a 

Consequently  part  (2)  equals: 


Evaluating  part  (3) 

This  part  contains  integrals  of  the  form: 


< 


I T2  exp(-  T2)  exp(-  ft)dt  . 
0 


As  T  «  e^(w  -  t)  rewriting  this  integral  one  gets: 

7T  ' f")  1  l2  '4  (T 

Let  a  “  — -r  .  As 
2e* 

jL  T  exp(-  (T  -  a)2)  -  exp(-  (T  -a)2)  -  2T(T  -  a)  exp(-  (T  -  a)2) 
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then 


e^w 


e*w 


I  T2  exp(-  (T  -  a)2)dT  m  J  f  «P("  (T  "  a>2>dT  ~  iT  exP("  <T  “  a>2> 


e*v 


C  W 

+  a  j  T  exp(-  (T  -  a)2)dT 


From  part  (1) 


e^w 


j  exp(-  (T  -  a) 2)dT  -  J^r(l  +  erf(e*w  -  a)j 


Evaluating 


e*w 

/ 


T  exp(-  (T  - 


a)  2)dT 


As 


^  exp(-  (T  -  a)2)  -  -  2(T  -  a)  exp(-  (T  -  a)2) 


then 


e^w 


/ 


T  exp(-  (T 


a) 2)dT 


-  |  exp(-  (T  -  a)  ) 


e^w 


eiw 


*/ 


exp(-  (T  -  ur)dT 


-  J  exp(-  (e^w  -  a)2) 


*sf(t 


+  er£(e^w  - 


a)) 


therefore 


Jt2  exp(-  T2)  exp(-  ft)dt  -  -y  exp^  -  fw^  |  (a2  +  i)  ^  (1  +  erf(e^w  -  a)) 
0  6 

-  {(e^w  ♦  a)  exp(-  (e^w  -  a)2) 


Thus  part  (3)  equals: 
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Evaluation  of  part  (4) 

The  integral  here  is  of  the  form: 


1 1  erf (T)  exp(-  ft)dt  . 
0 


As  T  »  e^(w  -  t)  changing  the  variable  of  integration  from  t  to  T  and  substituting 
for  the  factor  t  one  gets: 


e^w 


7  / 


erf (T)  exp(-  ft)dT  +  ~  j  T  erf(T)  exp(-  ft)dT  . 

e 


A 


J  V. 


B 


The  integral  in  component  A  has.  already  been  evaluated,  as 


00 

/ 


erf(T)  exp(-  ft)dt 


1 

J 


C  V 
/ 


erf(T)  exp(-  ft)dT 


which  from  part  CD  equals 

lerf(e^w)  -  expl 


1 

3 


n« 


(4 -^i)| 


then  denoting  this  by  part  (A)  then  component  A  equals  (w/neh  •  part  (A). 
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Rewritting  component  B  in  the  form 


y  exp(-  fw)  I  T  erf  (T)  exp^ 


dT 


as 


±  |t  .r«(T> 


erf  (T)  exp^yj  +  77  exp(-  T2)  exp^yj 


+  —r  T  erf (T)  e 
e* 


then 


— oc 

/ 


e^w 


T  erf  (T)  e  ttr 


_*  e^w 

T  erf (T)  exj/-^r)  +  J 

g4W  -CO 


\7) 


erf (T)  e 


3j}z 


e^w 


*-h  / 1  **><- 12> 


Evaluating  the  terms  in  the  above  expansion.  From  the  above 

e*w 

J  erf (T)  exp^y^dT  -  7  exp(£w)  •  part  (A)  . 


As 


e^w 


2 

77 


J  T  exp(-  T2)  exp^jdT  equals  exp^  “  fwj  j  T  exP(“  (T  “  “>2>dT 


where  a  -  f/2e^  ;  and  (from  part  (3))  as 


e*w 


J  T  exp(-  (T  -  a) 2)dT  equals  ^  (I  *  erf(eJw  -  a))  -  |  exp(-  (e*w  -  0)) 


L  _ 


X- 


then 


e^w 


2 

77 


M 


T  exp(-  T2)  expf-^T^  T 


yy  exp^|j  “  +  erf (e^a  -  w)) 

_ _ /_  /■-i„  _ 
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Consequently 


^-^■e^w  erf(e^w)  exp(fw)  +  ^  exp(fw)  •  part  (A) 

+  exp^£  ”  £wj(a/ir(  I  +  erf(e^w-a))  -  exp(-  (e^w-a)2^j 


and  therefore 


e 


T  erf (T)  exp(-  ft)dT 


17 


part  (A)  -  we^  erf(e^w) 


+  7^exp 


(fc  "  2£w){a/* 


a/*(\  +  erf(e*w-a)) 


l  2 

-  exp(-  (e2w-a)  ) 


f 


The  solution  of  part  (A)  is  therefore: 


part  (Aj) - jr  erf 


•  part  (A2)  - 


where 


and 


+  1JJ  6Xp(fe  "  2fW2)(Q7li  ('  +  erf  (*4w2  '  “))  “  exp(~  (\ 

part  (A,)  -  erf(eJw,)  -  exp^  -  fwj  ♦  erf^w,  -  , 

part  (A2)  -  erf^wj  -  exp^  -  fw^l  +  erf^Wj  - 


a 


n 


\ 
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Evaluating  part  (5) 

The  integral  of  interest  is  of  the  form: 


00 

/ 


t  exp(-  X)  exp(-  ft)dt 


As  T  «  e^(w  -  t)  this  can  be  written  as: 


CO  00 

w J  exp(-  T2)  exp(-  ft)dt  - y  j T  exp (-  T2)  exp(-  ft)dt 


B 


Utilizing  the  solution  given  in  part  (t)  then 


w  j  exp(-  T2)  exp(-  ft)dt  -  w  exp^|j-  -  fw^  -  erf^-£y  -  e^w^ 


To  evaluate  B,  complete  the  square  and  change  from  t  to  T  .  On  doing  this  B  becomes 


T  exp(-  (T  -  «)z)dT 


e^w 


A 


where  a  «  f/2e*  .  From  part  (3) 


/ 


.L  _  ~^2\  _ 

2 


T  exp(-  (T  -  a^JdT  «  J  exp(-  (e!w  -a)*)  -  (1  +  erf(e^w  -  a))  , 


e*w 


therefore 


J 


1  j T  exp(-  T2)  exp(-  ft)dt 


‘  ^  «*P^  "  fwj  |®*P<“  (®*W  ~  a> 


a/ii 


(1  ♦  erf(e*w  -  a)) 


3 
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Consequently  the  solution  of  part  (5)  is: 


w2  exp^  “  f"2)  X  °-5Jj  0  '  erf(°  "  ^2)) 

+  exp(4  “  fw2)  [exp(‘  (eiw2  “  a)2)  "  a/r(‘  +  erf  (eS  ' 
«,  e*p(4  -  fwi) *  °*5-J?  ('  '  erf(“  ‘  eS)) 

+  exp(S  ~  fwl)  {CXp("  "  °)2)  ‘  a/*('  +  erf  (elwl  " 


where  a 


f/2e*  . 


Evaluating  part  (6) 

The  integral  of  interest  is 


J tT2  exp(-  T2)  exp(-  ft)dt  . 
0 


As  T  -  e^(w  -  t)  this  can  be  written  as: 


■o  “ 

wJt2  exp(-  T2)  exp (-  ft)dt - y  J  T3  exp(-  T2)  exp(-  ft)dt 


■V 

B 


From  part  (3) 


wj  T2  exp(-  T2)  exp(-  ft)dt 
0 


V 


*)  {(°2  +  i) 


(1  +  erf(e^w 


-  l(e^w  ♦  a)  exp(-  (e^w 


-  w  •  Q 
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Appendix  C 


Rewritting  part  B  as 


i  exp^  -  fv)  / 


T3  exp(-  (T  -  a)2)dT 


where  a  -  f/2e*  and  considering  d/dT  (T2  exp(-  (T  -  a)  )  which  equals 
2T  exp(-  (T  -  a)2)  -  2T2(T  -  a)  exp(-  (T  -  a)2)  then 


exp(-  (T  -  a) 


_  —CO 

—CO  . 

2) dT  »  J  T  exp(-  (T  -  a)2)dT  +  a  J  T2  exp(-  (T  -  a)2)dT 


-  JT2  exp (-  (T  -  a)  ) 


From  part  (3) : 


T  exp(-  (T  -  a) 2) dT  =  i  exp( 


-  (e*w  -  a)2)  -  ~  (l  +  erf (e^w  -  a)) 


From  the  above 


J T2  exp(-  (T  -  a)  Z)dT  -  -  e*  exp^-  ^  -  fw^ 


therefore 


oo 

-j  J T3  exp(-  T2)  exp(-  ft)dt 


7  exp(fe  "  fw)  x 

X  exp(-  (e*w  -  a)2)  -  ^  (>  *  erfCe^  -  a)) 

-  ae*  exp^-  ^  -  fw^  •  Q  +  exp(-  (e*v  -  a) 


-  a)2)  -  (1  +  erf(e^w  _  a)) 
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Appendix  C 


Consequently  the  solution  to  part  (6)  is: 


and 


f 

a  *  — r  • 

On  * 


\ 
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Appendix  D 

THE  DISCRIMINATOR  CURVE 


A  quantitative  description  of  this  curve  is  required,  as  mentioned  in  section  4.2, 
so  that  the  tracking  bias,  ie  the  value  of  t  when  I(t)  -  0  ,  can  be  calculated.  This 
bias  is  used  to  correct  the  processed  height  values.  The  discriminator  curve  I(t)  is 
described  by  the  expression  I(t)  *  M(t)  -  aL(W  +  t)  where  M  and  L  are  the  middle  and 
late  gates  employed  by  the  tracking  gate,  t  is  the  time  offset,  ie  the  displacement  of 
the  centre  of  the  middle  gate  from  the  origin  of  the  waveform's  time  axis  while  W  is 
the  width  of  both  gates. 

Appendix  C  gives  the  quantitative  description  of  the  output  as  received  from  the 
DFB.  This  output  equals  the  convolution  of  the  instrument  transfer  function  with  the 
waveform.  As  the  instrument  transfer  function  It  has  a  'top  hat'  form,  ie 

t  .  '}  •T.'S,<T2 

'  oJ  ,<_V  T  51  t2 

and  the  middle  and  late  gates  are  also  'top  hat'  functions,  thus  M(t)  is  given  by  the 
solution  in  Appendix  C  where  Ti  “  T2  “  W/2  and  likewise  L(W  +  t)  is  given  by  the  same 
solution  where  Tj  *  W/2  and  t2  “  3W/2  .  By  using  the  above  expression  for  I(t)  , 
with  a  »  0.5  the  discriminator  curve  can  be  determined  given  W  . 

For  a  specified  pulse  width  and  value  of  a  the  tracking  bias  is  dependent  on  a 
(rms  waveheight) ,  X  (skewness) ,  E^  (attitude  error  angle)  and  W  .  Knowing  a^,  X 
and  Er>  W  is  selected  on  the  basis  of  minimising  this  bias  -  the  gate  widths  used  as  a 
function  of  a '  are  given  in  Table  2.  This  condition  ensures  that  the  tracking  gates 

always  sample  a  linear  portion  of  a  discriminator  curve  ensuring  high  sensitivity  and 
stability  of  the  ATU.  As  an  adaptive  tracking  gate  system  was  used  in  the  Seasat  1 
altimeter  the  three  gates  comprising  a  tracking  gate  are  of  equal  width  W  and  are 
contiguous.  By  employing  contiguous  gates  the  bias  variation  with  E^  is  minimised. 

Lastly  as  AGC  employed  in  the  altimeter  system  kept  the  integrated  DFB  output  of 
50  consecutive  waveforms  constant,  for  each  position  of  the  tracking  gate  during  the 
convolution  procedure  the  area  of  the  waveform  sampled  by  the  DFB  must  be  fixed  at  a 
constant  level.  This  area  determination  and  consequent  adjustment  were  carried  out 
numerically  in  the  simulation  process  described  in  section  4.3. 
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Table  1 

SEASAT  1  ALTIMETER  PARAMETERS 


Frequency 

13.5  GHz 

Pulse  repetition 
frequency 

1024  Hz 

Bandwidth 

320  MHz 

System  noise 
temperature 

3364  K 

Chirp  pulse 
width 

3. 2  ys 

Gain  control 

Automatic 

Compressed 
pulse  width 

3.125  ns 

Antenna  main 
lobe  gain 

40  dB 

Pulse 

compression  ratio 

1024 

Antenna 

polarization 

Linear 

Peak 

transmitted  power 

2  kW 

_ 

Antenna  diameter 

1  m 

Table  2 

GATE  WIDTHS  EMPLOYED  IN  CALCULATING  THE  TRACKING  BIAS 


Rms  waveheight 

(m) 

Gate  width 
(ns) 

5 

68.75 

4 

59.375 

3 

43.75 

2 

34.375 

1 

18.75 

0.5 

12.50 

45 


AGC 

ATU 


CW 

D 

dA 

DDG 

DDL 

DFB 


e 

EM 


E 

s 

E. 

js 

A 

e 

A 

e 

-z 

FM 

FT 

FT-1 


G 

GHz 


H 

H 

H. 

-l 

H 

-s 

h(x,y) 

HSWS 

I 

I(t) 

Ir(T) 

It(t) 

J 

k 


k. 

-l 

k 

-s 

K. 

J 

kHz 


L 

LO 

LSB 

M 

MHz 


LIST  OF  SYMBOLS  AND  ABBREVIATIONS 

automatic  gain  control 

adaptive  tracking  unit 

speed  of  light 

continuous  wave 

displacement  current 

element  of  area 

digital  delay  generator 

digital  delay  line 

digital  filter  bank 

eccentricity 

electromagnetic 

attitude  error  angle 

scattered  far  field  electric  vector 

contribution  by  the  jth  element  to  Eg 

unit  vector  of  the  incident  electric  field 

z  component  of  e 

frequency  modulated 

operation  of  the  Fourier  transform 

operation  of  the  inverse  Fourier  transform 

power  gain  of  the  altimeter  antenna 

gigahertz 

satellite  altitude 
magnetic  field  vector 

magnetic  field  vector  of  the  incident  EM  wave 
magnetic  field  vector  of  the  scattered  EM  wave 
surface  elevation  at  point  x,  y 
high  speed  wave  sampler 
unit  dyad 

tracking  law  output 
flat  sea  impulse  response 
instrument  transfer  function 
current  density 

magnitude  of  the  wave  number  vector 

wave  number  vector  of  the  incident  EM  wave 

wave  number  vector  of  the  scattered  EM  wave 

Gaussian  curvature  at  point  j 

kilohertz 

power  output  of  the  late  tracking  gate 

local  oscillator 

least  significant  bit 

power  output  of  the  middle  tracking  gate 

megahertz 
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LIST  OF  SYMBOLS  AND  ABBREVIATIONS  (concluded) 


MSB  most  significant  bit 

N  average  number  of  specular  points  per  unit  area 

n  unit  vector  normal  to  surface 


unit  vector  normal  to  surface  at  the  jth  element 

number  of  pulses  averaged 

nanosecond 

power  output  of  altimeter  transmitter 
power  received  from  altimeter  footprint  at  time  t 
probability  density  function 
slope  pdf 

joint  slope  height  pdf 
pulse  shape  distribution 

radial  distance  between  area  element  and  altimeter 
distance  measurement  in  plane  of  the  sea  surface 
principal  radius  of  curvature  at  point  i 
root  mean  square 

synchroniser,  acquisition  and  control  unit 
signal  to  noise  ratio 

time  measurement  with  origin  at  time  of  pulse  transmission 

width  of  middle  tracking  gate 

width  of  late  tracking  gate 

point  on  leading  edge  tracked  to 

tracking  gate  width 

waveheight  pdf 


number  of  DFB  gates  comprising  middle  and  late  tracking  gates 
A  increment  in  a  variable 

At'  height  resolution  size 

6(  )  Dirac  delta  function 

n  angle  of  incidence  at  each  point  j 

0  antenna  beamwidth 

X  skewness  coefficients 

X.  wavelength  of  impinging  radiation 

u  permeability  constant 

us  microsecond 

u _  moment  of  the  distribution  P(z;Ax) 

on 

o  backscatter  cross  section  for  a  unit  area  of  scattering  surface 

oR  rms  tracking  error 

Op  rms  width  of  the  distribution  describing  the  transmitted  pulse  shape 

ow  rms  waveheight 

o£  standard  deviation  of  surface  elevations 

o&%  standard  deviation  of  surface  slopes  in  the  x  direction 

t  time  measurement  with  origin  at  time  of  recaption  of  signal  from  mean  sea 

level 

v  angular  frequency 
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Fig  1  Schematic  of  the  Sea  tat  1  altimeter  tyitem 
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Fig  2  Geometry  illustrating  the  Interaction  of  a  spherical  pulse  with  a  flat  sea  surface 
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Fig  4 


llluitration  of  two  typical  scattering  surfaces. 

Top  -  smoothly  varying  surface  where 
Bottom  -  rough  surface  wbera  R  <  1 
(h(x,y)  -  surface  elevation  at  point  x,y) 
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Fig  6  Surface  exhibiting  multiple  scattering  and  shadowing 
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